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Abstract 

D 

We apply the nonlinear steepest descent method to a class of 3 x 3 Riemann-Hilbert problems introduced 
in connection with the Cauchy two-matrix random model. The general case of two equilibrium measures 
supported on an arbitrary number of intervals is considered. In this case, wc solve the Riemann-Hilbert 
KH problem for the outer parametrix in terms of sections of a spinorial line bundle on a three-sheeted 

Ricmann surface of arbitrary genus and establish strong asymptotic results for the Cauchy biorthogonal 
polynomials. 
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1 Introduction 

In this paper, we study asymptotic behavior of a class of 3 x 3 Riemann-Hilbert problems motivated 
by the recently introduced two-matrix random model [7] . The model consists of two random Hermitean 
positive-definite matrices Mi , Mi of size n x n with the probability measure 

&H{Mt,M 2 ) = d f' dM ', e -NTr m M l))e -NTr { V(M 2)) ^ 

PV ' ; Z n det(Mi + M 2 ) n K ' 

where U, V are scalar functions defined on K + . The model was termed the Cauchy matrix model because 
of the shape of the coupling term. Similarly to the case of the Hermitean one-matrix models for which 
the spectral statistics is expressible in terms of appropriate biorthogonal polynomials [20], this two- 
matrix model is solvable with the help of a new family of biorthogonal polynomials named the Cauchy 
biorthogonal polynomials [6]. 

The Cauchy biorthogonal polynomials are two sequences of monic polynomials (p J (a;))°° =0 , (qj(y))jt 
with deg pj =deg qj = j that satisfy 

e -N(U{x)+V{y)) 

Pj(x)q k {y) dxdy = h k 5 jk , Vj, k > , h k > . (1-2) 

xR + x + V 

These polynomials were studied in [9] in relation with the spectral theory of the cubic string. 

In yet another application, in analogy to the moment problem approach to the Camassa-Holm peakons 
[3, 4], the predecessors of the Cauchy biorthogonal polynomials were used to study the peakon solutions 
of the Degasperis-Procesi wave equation [19, 18]. More generally, the Cauchy biorthogonal polynomials 
are expected to play a role in a variety of inverse problems for the third order differential operators [21]. 

The main features of Cauchy BOPs can summarized as follows: 

1. they solve a four-term recurrence relation; 

2. their zeroes are positive and simple; 

3. their zeroes have the interlacing property; 

4. they satisfy Christoffel-Darboux identities; 

5. they can be characterized by a pair of 3 x 3 Riemann-Hilbert problems; 
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6. the solution to the corresponding Riemann-Hilbert problem yields all kernels for the correlation 
functions of the Cauchy two-matrix model; 

Items 1-5 above have been addressed in [6] while item 6 was explained in [7]. In the present paper, 
we apply the Deift-Zhou steepest descent method to the asymptotic analysis of the Riemann-Hilbert 
problem with the view towards applications to the biorthogonal polynomials and the spectral statistics 
of the Cauchy two-matrix model. 

The paper is organized as follows. In Section 2, we set up a Riemann-Hilbert problem that char- 
acterizes the biorthogonal polynomials and that is essential in evaluating the correlation kernels for the 
associated matrix model. 

In Section 3, we recall the results of [2, 7] where a relevant potential theory problem was set up and 
solved. The resulting equilibrium measures are the key ingredient in the construction of the g-functions 
that pave the way for the Deift-Zhou steepest descent method. 

The central section of the paper is Section 4, which deals with the nonlinear steepest descent analysis. 
We consider the general case in which the two equilibrium measures are supported on an arbitrary number 
of intervals. This prompts the use of a higher genus three-sheeted Riemann surface. Note that a Riemann 
surface of a similar structure was recently used in [22] to study Hermite-Pade approzimations of pairs of 
functions that form generalized Nikishin systems. 

We solve the Riemann-Hilbert problem for the outer parametrix in terms of sections of a spinorial 
line bundle in the spirit of [8, 5]). Much of the effort goes towards showing that the model problem for 
the outer parametrix always admits a solution (Proposition 4.4, Theorem 4.1). 

It is perhaps worth mentioning that in [13] (Section 8) the authors approached a similar problem of 
solving a 4 x 4 RHP for multiple orthogonal polynomials that arise in the analysis of the two-matrix 
models with interaction e -ArTr ( MlM2 ) instead of the Cauchy interaction in (1-1). Though similar, their 
approach differs in that they use the theorem on existence of meromorphic differentials on a Riemann 
surface without providing explicit formulae in terms of Theta functions, contrary to the present work. 
Also we work with arbitrary two measures without being restricted to the choice of a quartic potential. 
Section 5 uses the asymptotic analysis in the previous sections to discuss universality results for the 
spectral statistics of the two-matrix model. We find that individual spectral statistics exhibits the same 
universality phenomena as in the one-matrix model. We expect that the Cauchy two-matrix model might 
produce new universality classes in the case when supports of both equilibrium contain the origin. This 
case, however, is not considered in the present paper due to assumptions on the potentials (Assumption 
1). Relaxing this assumption requires a generalization of the potential theory problem studied in [2]. 

The Appendices contain a discussion of a specific genus zero example (Appendix A) as well as notations 
and essential information regarding Theta functions (Appendix B). 

Acknowledgements. This paper was completed at the Banff International Research Station. We thank 
BIRS for the hospitality and for providing excellent work conditions. We also thank Dima Korotkin for 
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helpful discussions concerning spin bundles. M. B. and J.S. acknowledge a support by Natural Sciences 
and Engineering Research Council of Canada. M.G. is partially supported by the NSF grant DMS- 
1101462. 



2 Riemann-Hilbert problem for Cauchy biorthogonal polynomi- 
als 

In the case of ordinary orthogonal polynomials, a well-known characterization in terms of a Riemann- 
Hilbert problem was obtained in [15]. Here we present a similar characterization following [6]. 
For symmetry reasons, we define 



V^z):=U{z), V 2 (z):=V(-z). 
Following [2] we consider potentials Vj that are subject to 
Assumption 1 The potentials Vj(z) satisfy: 



Vj(z) is a real analytic function on (- 

the growth- conditions 



Y7+1T 



Vj(x) = -ctj In \x\ + 0(1) as x -> 0, a 3 > 1, 



lim 



x->-(-iy+ 1 oo In \x\ 

• the derivatives Vj(z) are meromorphic on a strip containing the whole real axis. 
Example 2.1 The typical examples are potentials of the form 

Vj(x) — —a-j In |i| + Pj(x) , a,j > 1 
where Pj are real polynomials such that lim a ._ > /_\j+i 0O Pj(x) — +oo. 

The relevant Riemann-Hilbert problem that characterizes CBOPs is 
Problem 2.1 Find a matrix T(z) such that 

1. r(z) is analytic in C \ R, 

2. T(z) satisfies the jump conditions 



T(z)+=T(z). 



1 e ~ NV ^ 
1 
1 



r(4 = r(4 

where the negative axis is oriented towards — oo, 



1 

1 e~ NV ^ 
1 



z e 



z e 



(2-1) 



(2-2) 



(2-3) 



(2-4) 
(2-5) 



3. at z = oo 

T(z) = [ I - O 

4- near z = 



z n 
1 
i 



(2-6) 



r(z)=[0(l),0(ln|z|),0(ln 2 |z|)]. (2-7) 

Remark 2.1 TTie growth condition at z = m egs. fS-7j can oe replaced by 0(1) if the densities e~ NVj 
vanish at x — 0. In /act, a/ier Assumption 1 is put in place, we have e~ NVj ^ = O (|a;| ajAr ) and this 
(using Plemelj formula for the local model ofT at z = 0) implies T(0) = O(l). 

In [6] (Section 6.2) it was shown that (adapting to the present notation x ) 
I» = (2-9) 
MZ) 2^7 R+ ~z W?k + L (y-z)(x-y) 

2l7rPn - liz) 1 + x- Z JJ R _ 2i,(y-z)(x-y) + ~toT 

(-) n (2m) 2 , , {-) n 2in f Pn Jx)e- v i^dx (-)" f f p^Ax^e^^^-^^dxdy 

- Pn -l{Z) 



h n -i K-x J R+ x-z h n -i J R+ J R (y-z)(x-y) 

where p n are certain polynomials of exact degree n described in [6], Wp*(z) = J R c y-l^ an d au 
integrals are oriented integrals. We will study the asymptotic behavior of the solution of Problem 2.1 in 
various regions of the complex plane for 

NBN^oo, n:=N + r,reZ, (2-10) 

where the integer r is bounded. 

3 The q functions 

An asymptotic treatment for the RHPs for T along the lines of the nonlinear steepest descent method [10], 
requires that we normalize the problem with the use of an auxiliary matrix constructed from equilibrium 
measures that minimize a functional described below (see [2], [7]). 



In loc. cit. the matrix that we denote here by T corresponds to T; the precise relation is 

T(z) =diag(l,-27ri,-(27r0 2 )f( z )diagfl,--^,-— *— ) (2-8) 

\ 2m (2m)- 1 J 
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Consider the space of pairs of probability measures p,\ , p, 2 supported on R + ,R_ respectively. On this 
space we define the functional 

S , [d ( ui,d^ 2 ] := / Vi(x)d/J,i(x) + f [ dp, 1 (x)dp 1 (x') In - — — - + 

JR+ JR + JM + \X — X \ 

V 2 (y)dfi 2 (y) + I I dp 2 {y)dp 2 {y')\n ' 



\y - y'\ 

d/ii(x)d^ 2 (y)ln|x-y| (3-1) 

The minimization of such a functional was studied in a more general setting in [2], it is related to a 
similar problem for a vector of measures of Nikishin type [24] . 

Theorem 3.1 (sec Theorem 3.2 in [2]) 

Under the Assumptions 1 there exists a unique pair of densities pi,p 2 that minimizes the functional 
(3-1). Moreover (see Fig. 1 ) the supports consist of a finite union of compact intervals and the densities 
Pj are smooth on the respective supports; 

Li L2 

supp(pi) = |J At C K+ , su PP (p 2 ) = |J Bi C R_ . (3-2) 



Remark 3.1 It has been recently proven by one of the authors and A. Kuijlaars that the growth condition 
near the origin in (2-2) can be simply disposed of. The properties in Theorem 3.1 are still valid except that 
the support of equilibrium measures may contain x — 0, in which case Pj(x) — 0(a;~t). This behavior is 
crucial in deriving a new type of universality near x = and will be part of a forthcoming publication. 

Theorem 3.2 (Theorem 5.1 in [7], see also [2]) 
The shifted resolvents 

r« : = - Wl + , ye> : = w 2 - YL±m t where Wj{z) := f EM dx , 

3 3 J x — z 

y(D + y(2) + y(0) =Q _ (3 _ 3) 

are the three branches of the same cubic equation in the form 

E(y,z):=y 3 -R(z)y-D(z)=0 , (3-4) 
where R(z), D(z) are certain functions analytic in the common domain of V{ and V 2 . 
As a corollary of Theorem 3.2 we deduce the 
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Figure 2: The jumps of the g^ and g^ 2 -* functions in the gaps. 



Corollary 3.1 For generic real analytic potentials the densities of the two equilibrium measures p\,p2 
vanish like a square root at the endpoints of each interval Ae,Bi in the support of the spectrum. 

Proof. It follows immediately from Cardano formulae, since the densities are related by Theorem 3.2 
to the jump-discontinuity of the branches of the pseudo-algebraic curve (3-4) which have in general 
square-root type singularities corresponding to simple zeroes of the discriminant A = 4R 3 (z) — 27 D 2 (z). 
Q.E.D. 



Definition 3.1 (g— functions) The g-functions are defined as 



(1) W 



(x) ln(> - x)dx ; g (2) (z) := / p 2 (x) ln(z - x)dx ; g (0) (z) + g^{z) + g^(z) = (3-5) 



where g^ is defined as an analytic function in the domain T>\ := C \ [ao,oo) and g^ is analytic in 
T>2 := C \ (— oo, bo] while g^ is analytic in T>o '■= *D\ H T>2 = C \ ((-co, bo] U [ao, oo)) 

Remark 3.2 We remind the reader that the integral defining g( 2 ' is an oriented integral. The orientation 
of the negative half-axis that we use may seem unusual. However, it is consistent with the one used in 
studies of multiple orthogonal polynomials on a collection of radial rays in the complex plane. 



Definition 3.2 The right and left cumulative filling fractions are defined as: 

C a 2i-1 



Pi{w)Aw, I = 1, 



p 2 (w)dw, 



(3-6) 
(3-7) 



b 



Note that e Ll = 1 while cr i2 = — 1. 

The variational inequal- 
ities that characterize the 
equilibrium measures pj [23], 
together with the Definition 
3.1 translate into the follow- 
ing theorem. 



B 3 B 2 B 2 Bi B i 


Ai 


A! 


A-2 


bi 6 


a B 


0.1 





Figure 1: The supports of the equilibrium measures pi,p 2 
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Theorem 3.3 (Analyticity properties for the g functions) 

The following properties hold 

1. gW is analytic in T>\ and has the asymptotic behavior gW = Inz + 0(z^ 1 ), z — > oo; 

2. g*- 2 -* is analytic in T>2 an d has the asymptotic behavior g' 2 ) = — lnz + 0(z^ 1 ), z — > oo; 

3. q(°^ is analytic in T>q and has the asymptotic behavior = 0{z~ r ), z — > oo; 

4- on the right cuts Ae the function 3 (g^\x) — g^\x)j is decreasing and there exists a constant 
7 + such that 

Q i i\x)-g^(x) + V 1 (x)+ 1+ = (3-8) 

5. on the right gaps Ae : 

Q+Hx) — g^^a;) = ~2iitee = constant £ iR, (3-9) 
(gV 0) (x) - g W (x) + V x (x) + 7+) > 0, (3-10) 

6. on the left cuts Be the function 3 (a^\x) — Q^\x)j is decreasing and there exists a constant 7_ 
such that 

gi 2) (x) - gf(x) + V 2 (x) + 7 _ = 0; (3-11) 



7. on the left gaps Be 



8. on the gaps Be [J At 



(2) (2) 

Q , (x) — Q_ (x) = —2inae = constant £ zR, (3-12) 
( f (x) - ^(x) + V 2 (x) +7-) > 0, (3-13) 



Q f(x)-^(x) = \ X&A J (3-14) 

Imoe x £ Be 



The proof is based on the facts that pj are positive densities, the variational (in)equalitics that arise from 
the minimization of the functional (3-1) [23, 2] and elementary complex function theory. We leave the 
details to the reader. 

Definition 3.3 If the statement that functions in items 4- and 6. above are decreasing is replaced by the 
requirement that their derivatives are strictly negative the inequalities in (3-10, 3-13) are strict then the 
potentials are said to be regular. 

Because of the expressions (3-10, 3-13) we make the following 
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Definition 3.4 Define the effective complex potentials by 

Vl {z) := V^z) - (1) W + (O) W +7+, (3-15) 
<p,(z) := V 2 {z) - +g( 2 )(z) + 7 _. (3-16) 

The form of the g-functions (Definition 3.1) implies some important inequalities and equalities for the 
effective potentials that we presently explore. 

Theorem 3.4 The effective complex potentials tp l , tp 2 satisfy the following properties 

1. For any £ there exists an open neighborhood of (a,2t-2, ci2e-i) for which the real part ofip 1 is negative 
away from the cut At; 

2. For any I there exists an open neighborhood of (62^-2, for which the real part of <p 2 is negative 
away from the cut Bg . 

Proof. It suffices to consider any of the right cuts, say, At- Then because of the identity 0^' = and 

(2) (2*) 

from the properties specified in Theorem 3.3 (since g\ — 01 =0) we have 

0i i) -0 w =-0f+0 w => 0i i) +0r=0 w +0 (o) - (3-17) 

More importantly, we see from (3-8) above that 

(1) (1) (1) (0) to 1 CA 

0+ - 0- = 0+ -Q+ ~ Vl - 7+ = ~<Pi+ = (3-18) 



= - [£' - W -V ± - 7+J = <Pi- (3-19) 

Equation (3-18) implies that the jump AgW is the left boundary value of the analytic function —<pi, while 
equation (3-19) represents the same quantity as the right boundary value of the analytic function cp±. 
The condition of strict decrease appearing in Theorem 3.3 then implies in view of the Cauchy-Riemann 
equations that Sfyjj strictly decreases as one moves perpendicularly away from the cut. Q.E.D. 

Since the three g-functions are antiderivatives of the resolvents Wj (3-3) we deduce that there must 
be two constants 71,2 such that 

(1) W = 7! + - f F«(C)dC ; 0< 2 >(z) = 72 - - f r( 2 >(C)d C , (3-20) 

6 Ja 6 Jb„ 

where the respective integrals are performed within the simply connected domains specified above and 
the constants 71,2 are chosen so that asymptotically 

g (1) (z) =\nz + 0{z- 1 ), q (2) (z) = - In z + Oiz- 1 ). (3-21) 

It thus appears that the three g-functions are intimately related to the three branches of the integral 
J ydz over the three-sheeted Riemann surface (3-4). 
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4 Deift— Zhou steepest descent analysis 



We will follow a well-established approach to asymptotic analysis of Riemann-Hilbert problems pioneered 
by Deift and Zhou in [11] and referred to as the nonlinear steepest descent method. 

4.1 Modifications of the Riemann-Hilbert Problem 

We proceed to introduce the contours indicated in Fig. 3: the disks D a around each endpoint a € 
{ctj, &j}, = x,.., of the cuts have sufficiently small radii so as not to include any other end points. The exact 
position and shape of the upper and lower arcs joining two end points are largely irrelevant; they should 
lie in the region where the real parts of the effective potentials (Definition 3.4) are negative. We will 
name the part outside all disks and all lenses the outer region. 

Let us define a new piecewise analytic matrix function (only nonzero entries are indicated) 



M_ 



6©& 




Figure 3: The final steps in the modification of the RHP. 



( 



1 e -"V! 



1 



1 



) 




Figure 4: The modified jumps for r o on a cut of the right side of the spectrum. 



10 



Figure 5: The modified jumps for T on a cut of the left side of the spectrum. 



r„ := < 



1 

3 w,W 1 



1 



1 

3 NV 2 (z) 1 



1 

_ e NV 2 (z) 1 



in the outer region 
in the upper half of each lens on the right cuts 
in the lower half of each lens on the right cuts 
in the upper half of each lens on the left cuts 
in the lower half of each lens on the left cuts 



(4-1) 



As a consequence the jumps of the Riemann-Hilbert problem are modified as indicated by the jump- 
matrices in Figs. 4, 5. We now define 

T(z) := C- 1 T Q {z)G(z)C 1 , G{z) := diagte-^ 1 ' « e - N ° m ^, e'^W) (4-2) 

(4-3) 



C-y ' 



e 3 



a 7V- 



-N 



The matrix T will be the main object of our interest from this point on. Recall that gy' = (— l) 3 lnz + 
0(z _1 ), for j = 1, 2, and that n = N + r. Then T satisfies a new Riemann-Hilbert problem: 



T(z) = (l + O(z- 1 ))z diag( ^ '- r) z^oo, T+(z) = r_(*)M(«) 



(4-4) 



where M(z) takes different forms depending on the arc considered. Note that T is bounded at z = 
because of Assumption 1 and Remark 2.1. Recalling the definition of the effective potentials (Def. 3.4) 
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-1 



1 




-1 




1 





+ 



+ 



Figure 6: The RHP for the outer parametrix 



we have (using the notation A/ = / + — /_, Sf = f + + /_) 

1 



M = M { ± " 9ht) 



M = M. 



(left) 



1 



1 



, on the upper/lower rims of a right lens, 
on the upper/lower rims of a left lens, 



(4-5) 



M = M, 



(right) 




_ Jv(v 1 - fl <! ) + fl y o, +7 + ) 



-1 



on the right cuts, 



M = M, 



Aleft) 







on the left cuts, 



M = M {r%9ht) = 

gap 



on the right gaps, 



M = M {left) = 

gap 



-2HvN<r e N -fl™ - V 2 —f. ) 



(4-6) 



(4-7) 



(4-8) 



(4-9) 



on the left gaps. 

Observe that the off-diagonal entries of (4-5), (4-9), and (4-8) are exponentially small as N — > oo 
because of the signs of the real parts; in particular all their L p norms (except forp = oo) are exponentially 
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Figure 7: The contours of the residual RHP for the error term and the orders of the jumps matrices. 

small. Yet, near the endpoints Oj, bj of the support of the equilibrium measures, the off-diagonal entries 
tend to 1. However, the new RHP that effectively amounts to setting them to zero is a key ingredient in 
the final approximation we need. It is, in fact, the main new ingredient of the paper from the point of 
view of the nonlinear steepest descent method. More explicitly, we formulate 

Problem 4.1 (Outer parametrix) Find a 3x3 matrix ^(z), analytic inT>o := C\((— oo,&o] U [ a 07°°)) 
and with the following properties 

(i) the jumps indicated in Figure 6; 

(ii) the growth conditions at z = oo and near an endpoint z = a are, respectively 



¥(z) = (l + o(-\\l 1 , *(*) = C ((*-«)"*) , «6Kilk.. (4-10) 



The final approximation involves the solution of the exact RHP for T within disks around the endpoints 
that we have indicated above (see Figure 3). This is called the local RHP and is formulated below 

Problem 4.2 (Local parametrix) Let D = D a be a disk (previously introduced) around any of the 
endpoints a of the supports of pi, p%. Find a piecewise analytic 3x3 matrix V{z) = V a (z) such that 

(i) V{z) is bounded for z 6 D uniformly with respect to N ; 

(ii) withinD, V(z) satisfies the jump- conditions for T with jump-matrices (4-9, 4-8, 4-1, 4-5); 
(Hi) on the boundary dD, V{z)^ [z) = 1 + o(l), N — > oo, uniformly in z G <9B. 

Suppose now that we have managed to solve Problems 4.1, 4.2 and let us define the matrix 

Then the "error" matrix 8(z) := T(z)T (z) solves the RHP 

£ + = £_(l + G), £(z) = l + 0(z- 1 ) , (4-12) 
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where the jumps are supported on the boundaries c?O a , on the rims of the lenses and gaps outside the 
local disks. A direct standard inspection ( see, e.g. [10]) shows that G tends to zero as N — > oo in all 
LP norms p £ [1, oo] and hence £(z) is close to the identity matrix 1 uniformly on P 1 . If the equilibrium 
problem is regular (in the sense of Remark 3.3) then the bounds on the jumps of the error term £ (and 
hence on the error itself) are as depicted in Fig. 7. 

4.2 Outer parametrix 

The construction of the solution to the RHP 4.1 uses Theta- Functions associated to an algebraic curve 
studied below. Solutions to problems of this kind can be derived in several ways [17, 8]. While we will 
eventually write down "explicit" formula? involving Theta functions, these formulae are not numerically 
effective unless the underlying Riemann surface has genus zero. However, we do need to prove existence 
of a solution. Thus, our strategy will be to produce theta-functional expressions for the solution to 
Problem 4.1 and subsequently use results on bordered Riemann surfaces from [14] to ensure solvability 
of the RHP in terms of the proposed expressions. As an example, we will also provide the corresponding 
(explicit) expressions for the case of genus zero. 

4.2.1 The abstract Riemann surface 

The title of this section refers not to the Riemann surface defined by the pseudo-algebraic equation 
E{y,z) = (3-4), but to an abstract Riemann surface C described below. 

Let L\ , L/2 be the number of intervals of the supports of the measures p\ , pi ; note that endpoints of 
cuts correspond to zeroes of odd multiplicity of the discriminant equation A(z) := 4i? 3 (z) — 27D 2 (z) , 
and the assumption of regularity of the potentials translates into the requirement that all zeroes of odd 
multiplicity of A(z) are real and simple. 

We consider now the abstract Riemann surface obtained by gluing together three Riemann spheres 
parametrized by the variable z; the Riemann sphere labelled 1 is slit along the support of p\ and glued 
there with the middle Riemann sphere labelled 0. The latter is subsequently slit also along the support 
of pi and glued across it with the Riemann sphere labelled 2, cut along the support of pi. 

The resulting Riemann surface £ is a compact surface of genus 

g = genus(C) = L x + L 2 - 2 (4-13) 

as follows from the Riemann-Hurwitz formula. The endpoints of the intervals Ae, £ = 1, ...£1 and 
Bi, t = 1, ...,L% are branch points of order 2, i.e. the local coordinate is \J z — xq, with xq one of the 
endpoints. The local coordinate around the three points at infinity 001,000,002 is 1/z. 

This point of view allows one to think of y^ 1 * 2 ) appearing in Theorem 3.2 as the three branches 
of a locally analytic function y over the Riemann surface C; indeed near a branch point xq, due to the 
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Figure 8: The Hurwitz diagram of the abstract spectral curve L: the vertical dotted lines represent the 
identification of points. 

assumption about A and due to Cardano formulae, the function y has a square-root Puiseux expansion 
in z — xq, and thus is a well defined function on the curve C. 

If V{, V2 are meromorphic functions then y is a meromorphic function on £, otherwise y will in general 
have other isolated singularities or will be defined only in a strip around the three-copies of the real axis 
in L. 

We will use the basis of the canonical homology of C indicated in Fig. 9. 




Figure 9: Our choice of the canonical homology basis on the abstract Riemann surface C; the solid lines 
represent arcs on Sheet 2, the dotted ones are arcs on Sheet 0, and the dashed ones are arcs on Sheet 1. 
In the example there are 5 total cuts and hence the genus is 3. We can declare that the curves lying on 
the same sheet are the a-cycles whereas the curves lying on two different sheets are the corresponding 
/3-cycles. 

Proposition 4.1 The curve C possesses a natural antiholomorphic involution defined as the complex 
conjugation of each of the three sheets. The finite gaps provide g — genus(C) point-wise invariant non- 
trivial cycles. Thus C can be realized as a double of the bordered Riemann surface obtained by identifying 
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Figure 10: Depiction of the lifted contours j x , j x i , 7 a , j y / . Their X-projection is just a closed path (in 
the picture we have chosen circles) intersecting the real axis only at x and one of &o or ao. The points 
x' ,y' belong to the gaps and by definition the curves "/ x >,"/ y ' he on one sheet and are actually homologic 
to the a-cycles chosen before (see Fig. 9 and its caption). 

the upper half of sheet 2 with the lower half of sheet along the left cuts, and the lower half of sheet 
with the upper half of sheet 1 along the right cuts (see Chapter 6 in [14])- 

Proposition 4.1 shows that the properties of C are very similar to the properties of a hyperelliptic curve 
defined as w 2 = P(z) for a real polynomial P of degree 2g + 2. This fact will be used to obtain information 
on the Theta divisor of the Jacobian of C when constructing the outer parametrix. 

4.2.2 Multiplier system \ an d spinors 

Recall the definition of domains 2?o,i.2 (Def. 3.1). By lifting points in the complex plane z one can 
construct three sections pj : T>j — > C such that: p\ is analytic in T>\ = C \ UAe, Pi is analytic in 
T>2 = C \ UBg and po is a analytic in T>q = T>\ n T>2 with the following boundary values 

P2{x)± = Po {x) T , xeB:= uf^Be (4-14) 

Definition 4.1 Let x £ M.+ (x £ the contour j x : S 1 — > £ is defined as the unique lift to C of the 

closed path on C that contains the leftmost (rightmost) point ao € A (bo G B) and intersecting the real 
axis only at ao and x (see Fig 10). The lift is accomplished by using the map p\ : T>\ — > C for the part in 
the upper half-plane (or pi : T>% — > C if x £ M.— ). The lower part of the path is lifted using po if x belongs 
to a cut or pi (p 2 ) when x belongs to a gap. 

By this definition j x defines a closed cycle in the homology of C. Moreover, ^ x span the whole 
homology as x ranges through R. 
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We define the following vectors (characteristics) in the Jacobian of the curve C 



A = 2inN 



°~L 2 -1 



(4-15) 



where the cumulative filling fraction have been introduced in (3-6, 3-7). Note the linear dependence of 
the vector A on TV . 

Given any two vectors A = (A,-),B = (Bj) € C 9 one can define a character, namely, a homomor- 
phism x '■ 7r i('C) ~~ ^ C* (with 7Ti the fundamental group of the Riemann surface) by extending the values 
it takes on the basis of the homology, xi a j) = &Aj i xifij) — eBj ■ For A, B in (4-15) we have 



x{otj) 



e 

N2iir<7 



1. 



(4-16) 



A meromorphic spinor ?/> for the character x is a multivalued half-differential that acquires the 
factor x{l) under analytic continuation along a closed contour 7. Any half-differential </> at a point p can 
be expressed in a local coordinate C as <^(p) = /(OV^C- I n particular, since z is a local coordinate away 
from the branch points, we can write vp(p) = f(z)\/dz, where p € C is a preimage of z. Recalling the 
definition of the three sections Pj(z) given above we have 



ip(pi{x))± = x(lx)^(po(x)) T , x e A , 
ip(P2(x))± = x(lx)ip(po(x)) T , x e B , 

ip{po{x))+ = x(7*rV(po(20)- . x e IF 



ip(pi(x))+ = x (j x )ip{pi(x))- , 
i>(P2(x))+ = x(7»)^(P2(*))- . (4-17) 
+ \(AUB) (4-18) 



In addition, we will use the spinor vdz defined on C slit along the top of the cuts A in sheet 1 and 
bottom on sheet 0, and along the top of B in sheet 2 and bottom in sheet 0. By definition \J dz satisfies 
the boundary conditions 

(4-19) 
(4-20) 
(4-21) 
(4-22) 



dz(px(x))± = ±Vdz(p Q (x)) T , , x € A 
/ dz(p 1 (x)) + = Vdz(p-i(x))- , 

xeB 



/dz(p 2 (x))± = ±Vdz(p (x)) 
f&z{jp 2 {x)) + = Vdz(p 2 (x))- , xeR\B 



Note that the jump relations (4-19, 4-21) imply that Vdz is globally defined on a double cover of C 
branched at the ramification points (to be distinguished from the branch points) since a 27r-loop around 
one such point (i.e. a 47r-loop around x — c for c a branch point) yields a transformation \fdz — > — vdz. 
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If we define a multivalued function f(p) := it follows from the above that its composition with 

y/dX(p) 

the three sections pj (Sec. 4.2.2) satisfies the following boundary-value conditions 

f(Pi(x))± = ±x(7x)/(Po(aO)=F > x e A 

/(pi(a;)) + = X (l x )f(pi(x))- , xeR+\A (4-23) 
f(P2(x))± = Tx(lx)f(po(x)) T , x G B 

f(p 2 (x))+ = x(7*)/(P2(a:))- , xeM_\6. (4-24) 

For our choice (4-16) of x the above relations take the form: 

f{pi(x))+ = f( Po {x))- , f(po(x))+ = -f(pi(x))- , xeA 

f(pi(x))+ = e 2 ^7( Pl (x))_ , x G (a 2j -i,a 2j ) (4-25) 

/(p 2 (a;)) + = -f(p (x))- , f( Po {x))+ = f(p 2 (x))- , xeB 

f( P2 (x))+ = e 2M 'f(p 2 (x))_ , x G (bu-ubu). (4-26) 

Moreover, near a branch point x — a, since the local coordinate on the curve C is \J z — a, the functions 
fj(z) :— f( P j(z)) behave as 

f 3 (z)=0((z-a)-i). (4-27) 

All the above discussion amounts to the following 

Proposition 4.2 Let ip be a meromorphic spinor for the character x and let fAC) := ^ p AC)) = ^ 27j eri 
the vector 

F(z):=[f 1 (z),f (z),f 2 (z)} , z€C\ (-00,60] U[a ,oo) (4-28) 

has the properties 

1 \ /^Ne 3 \ 

F(.r)_= F(.,-)_|-l 1 F(a:) + = F(a:)_ f e -^Ne 3 \ ; x G (<Z2j_i, G^j) 




1%-kN at 



, z g (b 2 e-i,b 2 e) 

(4-29) 



as lueZZ as F(a;) = C((x — a) *) for any endpoint a. 



This means that the jump conditions for the Problem 4.1 are already satisfied. The complete solution 
amounts to finding a spinor ip for each row such that it has an appropriate growth at the points above 
z = 00. Gcnerically, a spinor for a character x is uniquely determined, up to a multiplicative constant, 
by choosing a divisor of degree —1. This is a consequence of Riemann-Roch-Serre theorem [16, 1] for 
line-bundles (our line bundle is the tensor product of a flat line bundle defined by the character x and a 
spinor bundle, namely a line bundle whose square is the canonical bundle) We have 
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== | fr-liPl) fr-M fr-M | , /r := "4= , /r := -7= ( 4 " 31 ) 

Vaz Vdz 



Proposition 4.3 Consider the two sequences of spinors tpr^tpr, r £E Z for the same character \ (4~16) 
satisfying the following divisor properties (see (B.15) in the appendix) 

(ipr) >~(r+ l)ooi +roo 2 ; ($r) > -(r + l)ooi - oo + (r + l)oo 2 (4-30) 

Consider the matrix 

JriPl) JriPo) Jr(P2) 

fr-liPl) fr-l(Po) Jr-\(Vl 

Jr-l{Pl) /r-l(Po) /r-l(Pa) 

TTien 2J(^) solves a RHP with the jumps (4-29), the asymptotic behavior 

VO{z) =diag(C 1 ,C 2 ,C 3 )(l + 0(z- 1 ))diag(z'',l,z- r ) (4-32) 
and Problem 4-1 admits a solution if and only if C1C2C3 7^ 0. 

Proof. By Proposition 4.2 the jumps conditions are automatically fulfilled. As z — > 00, Pj(z) —> 00 » £ C 
by construction. Then 

f r ( Pl ) = KPz r (l + G(z- 1 )) J r {p )=Ki 0) l(l + O(z- 1 )) J r (p 2 )=Ki 2) z- r - 1 (l + 0(z- 1 )), 
^_i(pi) = ^V- 1 (l + 0(z- 1 )) ,/ r _ 1 (po)=%(l + 0(z- 1 )) ,/r-i(Pa) = ^ ) 1 2- r ~ 1 (l + 0(«- 1 )) 

This asymptotic behaviour follows from the divisor properties of iji r , ip r ; for example, since ip r has a pole 
at most of degree r + 1 at 00 1 where Vdz has a simple pole, / r must have at most a pole of order r at 
infinity. 

The last claim follows from the fact that the Problem 4.1 requires that the constants K r X \ K^2 1 and 
should not vanish. It is only in this case that the matrix 03 can be normalized by a left multiplication 
to behave like (1 -|-0(2:~ 1 ))diag(z r , 1, z~ r ). To prove the necessity of the last claim we have to show that 
if any of the constants Kr , K^\, vanish, then the Problem 4.1 is unsolvable. To see this we first 

point out that standard arguments show that if a solution to Problem 4.1 exists, then that solution is 
unique. Therefore the solution of Problem 4.1 exists if and only if the only row- vector solution F(z) of 
the RHP with the same jumps as 4.1 but with the asymptotic condition (4-10) at z = 00 replaced by 

F(z) = [0{z r -% Oiz-^M^ 1 )] (4-33) 

is the trivial solution F(z) = 0. If any of the constants above is zero, the corresponding row then provides 
a nontrivial solution precisely to this latter problem. This proves the necessity of the last claim. Q.E.D. 
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4.3 Theta— functional expressions 



The notation below is borrowed from [8] and the definitions can be found in [14] and are reviewed in 
Appendix B. 

The spinors tp r , ip r in Proposition 4.3 can be written "explicitly" in terms of Theta functions as follows 



VV-1 



6a(p~ 002) 
Q A (p - 001] 



r+1 



e 



(p + roo 2 - (r + 1)ooi)/ia(p) 



Oa(p- oo 2 ) 



-e 



(p + roo 2 — ooq — rooi)/iA(p)- 



(4-34) 
(4-35) 



6a(p- ooi) r Q A (p~ oo ) 

Here h A is a certain fixed holomorphic spinor defined in (B.10. We now recall that Qa(p~<i), as a function 
of a point p, has a simple zero at p = q and at other g—1 points whose positions are independent of q and 
depend solely on the choice of odd characteristic A. Also, by construction [14], the spinor Iia has zeroes 
precisely at the same g—1 points (and of the same multiplicity). It thus appears that t[ r has a zero of 
multiplicity r at 00 2 and a pole of order r + 1 at 00 1. The order of the pole at 00 1 may be smaller if the 
teimTAp) := 6 £ {p+ roo2 — (r+l)ooi) vanishes there. If this happens, the constant C\ in Proposition 
4.3 would be zero and the RHP unsolvable. Similarly, if 7^-1(002) = then ip r -\ has a zero of multiplicity 
higher than r, and consequently C3 = 0. Finally if T r —\(p) := O g (p + r inf ty 2 — 000 ~ rooi) vanishes 
at ooo, then the spinor "0r-i does not have a pole at 000 and thus C 2 in Proposition 4.3 vanishes. Note 
that 



(roo 2 — rooi) 



%(00i) = T r -l{c0 2 ) = %-l{00 Q ) = 6 

Therefore we have the proposition below 

Proposition 4.4 The Riemann-Hilbert problem 4-1 is solvable if and only if 



e 



(roo2 — rooi) 7^ 



(4-36) 



We now need to establish the nonvanishing property required by Proposition 4.4. We recall that in 
our case A £ iM. 9 and B = (4-15). We then can use the results on bordered Riemann surfaces contained 
in [14]. 



Theorem 4.1 For any A £ 2inR 9 , r £ Z, 

e 

Therefore, Problem J^.l is always solvable 



{roo 2 — rooi) 7^ 



(4-37) 



Remark 4.1 The theorem could be rephrased as stating that the monodromy data for the outer parametrix 
do not lie on the Malgrange divisor. In this case the Malgrange divisor is actually identifiable with the 
(0) divisor in the Jacobian of the curve C [17]. 
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Proof. Recall from Proposition 4.1 that L has an antiholomorphic involution and can be represented as 
the double of a bordered Riemann surface ([14], Chap. VI). 

The /3-cycles are homologous to cycles fixed by the antiholomorphic involution, since we can realize 
each of them as a path following the gap in R on one of the outer sheets and the same segment on the 
middle sheet. There are in total g + 1 cycles fixed by the involution; aside from the g /3-cycles, there is 
also the cycle which we denote by /3 s +i that covers the union of the segment [&o,ao] on the middle sheet 
with all the unbounded gaps in the three sheets (see Fig. 1). 

We now use Corollary 6.5 on page 114 of [14] which describes the intersection of the O-divisor with 
a (/-dimensional real torus in the Jacobian of the form * (with fj, = (/.if) 6 |Z ff a half-period). It 
consists of the image in the Jacobian of the set of divisors of degree zero that have a form X> ff -i — 
where 2?°_i is the divisor of degree g — 1, such that 22?°_ 1 is linearly equivalent to the canonical divisor 
and 2? s -i is any positive divisor of degree g — 1 invariant under the involution and such that £> g _i has 
(1 + 2fij) mod 2 points in (3j for j = 1, . . . , g. This means that if there existed a point A £ iR 9 such that 
9 q (0) = 0, then it would be the image of a positive divisor of degree g — 1 with exactly one point in 

each of the g gaps, a clear contradiction (see also Proposition 6.16 in [14]). This proves that e MM (o) 
cannot vanish for any value of A 6 zK 9 . 



To complete the proof, recall that 9 



(rcx>2 — rooi) is proportional to 6(A + roo 2 — rooi) . To 
prove that the latter expression is nonzero we only need to show that u(roo2 — rooi) £ iR g , where u is 
the Abel map. By our choice of a-cycles and normalization (B.l) , holomorphic differentials luj satisfy 
u>j(z) = —ujj(z). Therefore an image under the Abel map of any divisor of degree contained in /3 g +i is 
purely imaginary. Since 001,2 G Pg+ii the proof is complete. Q.E.D. 

The expressions in terms of Thcta functions for the columns of the outer parametrix ^ possess other 
interesting relations that were investigated in [5]. 

4.3.1 Genus case 

A particularly simple situation is the case in which the surface C is of genus 0, namely there are only two 
cuts, one in M + and one in R_ . In this case we can write quite explicit algebraic expressions for all the 
objects above (see Fig. 12 for an example). First, we introduce a uniformizing parameter t : CP 1 —> C 
in terms of which the meromorphic function z : L — > CP 1 is written as 

z:=u t + C+^ l + ^ l . (4-38) 

Here we have identified the three poles 000,1.2 with t = 00, 1, —1 respectively. Once the uniformization 
has been fixed, the expressions for the spinors above are extremely simple 



Vv (t + iy /dt _ (t + iy +1 

Jr ■ r~i r± \/ i . ' J r 



dz (t - l) r+1 V dz ' Jr (i - l)'V"o(< 2 - I) 2 - Mt + l) 2 -u 2 (t- l) 5 
? gr _ {t + i) r+1 [dt ? _ (t + l) r + 2 



(4-39) 



(t - l)' r +! V dz ' (t - l)V"o(< 2 - I) 2 - «i(i + I) 2 -u 2 (t- l) 2 
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The normalization is obtained by expanding near the three points at infinity 



fr 



i u 



r+i 



z r+l 2 r + 1 



t ~ 1 



t ~ 



fr- 



-r-l 



Z 



-l tto 



z r+l 2 r 
1 



(4-40) 



The solution to the RHP problem Problem 4.1 is 

fr(h(z)) 



t ~ 



Mt (z)) fr(t 2 (z)) 



*'(-) =<H^( (-",')' ^ '-'X ■»,,'« 2 r+ ') | fr-Mz)) f r -l(t Q (z)) fr-Mz))] • (4-41) 

Jr-l(h(z)) fr-l(to{z)) f r -l(t 2 {z)) j 

Here fi,2,3(^) are f-coordinates of sections pi, 2,3 discussed in Section 4.2.2. 
4.4 Local parametrix (solution of Problem 4.2) 

For the sake of completeness we spell out the form of the parametrix near a cndpoint of a cut where the 
density pi, 2 vanishes like a square root (see Corollary 3.1). 

4.4.1 The rank- two parametrix 

We need to solve the exact RHP in Fig. 11. This part is essentially identical to the established results 
in [12, 10]. We consider only the previously defined neighborhood D a of the right endpoint a := a 2 e-i 
of one of the intervals of p\. The modifications for the other cases are straightforward. Near x = a 
the effective potential tp 1 is a piecewise analytic function with the jump tp 1+ — ip 1 ^ — Ainei, where £ 
is the number of the gap that starts on the right of a and eg is the cumulative filling fraction (3-6). 



Let 1/3 be a function defined by 

M _ / (fi(z) + 2in(l-e t ) , 3z > 
nZ> \ cpx(z) + 2iTr(l + e e ) , 3z<0 

Then tp has an analytic expansion in 
Puiseux series of the form described be- 
low. 

Lemma 4.1 The locally analytic func- 
tion if has an expansion 

<p(x) = ^C(x-a)i(l + 0(x-a)) 

" C= hm ^M>0< 




Figure 11: The exact Riemann-Hilbert problem near a right 
endpoint in H. + , and in terms of the zooming local coordinate 



>a~ y/ a — X 

where the cut of the root extends on the 

left of a and the term 0(x — a) is analytic at a (i.e. has a convergent Taylor series) 
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Proof. Recalling the representation (3-20) for the g functions we have 



<p{z) 



-d£ 



(rW(o-y (0) (0) d £ 



(4-42) 



Now a small loop around a interchanges <H- Y^°\ which means that 



F (i)(£) _ y(o) (e) = ^~ a {2C + 0(£ - a)) (4-43) 
with the term 0(£ — a) analytic. The expression for C is obtained by recalling that by Def. 3.4 



^(z) = ^(z) = V((z) 



pi(s)ds 



p 2 (s)ds 



(4-44) 



Now pi(s) ~ Cy/a — s near s = a (with C > 0) due to (4-43) and <p' + — (p'_ = Amp\. The sign of C > 
is due to the fact that $t(p 1 = 3iip and it must be positive in the gap (x > a). Q.E.D. 



Remark 4.2 In a non-generic situation, one can only conclude that 

<p(z) = {z-a) k+ i{C + 0{z-a)) , C > 



(4-45) 



for some nonnegative integer k. This case is called nonregular and it is treated in [10]. The construction 
there applies identically in this case, but we prefer to concentrate on a detailed discussion for the most 
generic case. 



We define the zooming local parameter by the equation 

U*=N<p(z) 



(4-46) 



so that the fixed-size neighborhood of a in the z-plane is mapped conformally (one-to-one) to an homo- 
thctically expanding neighborhood of the origin in the £ plane with a diameter growing like N*. The 
choice of the root in (4-46) is such that the cut is mapped to M_ of £-plane. 

We then introduce the standard rank-two Airy parametrix [12, 10] as the piecewise defined 

matrix constructed in terms of the Airy function Ai(a;) as follows. If <x 3 = 



1 

-1 



define 



2/0 -Va\ e (§?2+^A^)<r 3 

y'o -y-J 



2/0 

y'o -j/2 



Kftt) := \/2^e-T J 



-Vx 

-y'x 

-2/2 2/i 

-2/2 y'i 

2/o 2/i 

2/o y'i 



-272 
-272 



\CS+iirNe,)tr 3 



§f2-i7rJVeJ<r 3 _ 



"3/1 

-3/i 

-2/2 2/i 

"3/2 2/i 



" J/2 |e^ ¥ 'i CT3 fc = 2 

"2/2 



eT^i 173 = fc = 3 



(4-47) 



2/o 2/i 
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where index k refers to regions in Fig. 11 and, for s = 0, 1, 2, 

y s :=w'Ai(w'0, uj = c 2 ™'\ 
Each ~RP k has the following uniform asymptotic behavior near £ = oo. 

<-±iirNe t )trati + 0(f~ 3 / 2 )). 



1 



(4-48) 



(4-49) 



where the ± depends on the half-plane (upper/lower) in which the asymptotics is considered. 

Thus, the final form of the local parametrix is defined as a matrix 1Z whose restrictions to the four 
regions in Fig. 11 are 



V 2 



1 -1 
1 1 



(4-50) 



Here 



(i) the prefactor F(z) solves a RHP on the left 



F{z). 



F(z). 



-1 

1 

F(z)+ = e~ 2i7rNe ^ 3 F(z)- , £(z) e M+ 



(4-51) 



(ii) 1Z satisfies the exact jump conditions of the RHP in Fig. 11 , except for the jump condition on the 
cut : 



-1 

1 o 



1 

-1 



(4-52) 



n+ = e -2i-KNe e cr 3n _ e 2vKNe e cr 3 £ £ R+ 

(iii) 7\L(£) = 1 + 0(N^ 1 ) uniformly on the boundary of the neighborhood D a . 
4.5 Rank-three parametrix 

Now we embed the rank-two parametrix from the previous section into a matrix 

n a {x) = n{x)®i=; 

Let define the parametrix V a (z) within the disk D a as 

V a (z) := 9(z)K a (z) 



TZ(x) 






1 



(4-53) 



(4-54) 



where 'J is the outer parametrix constructed earlier. The only point to address now is whether ^ (z)TZ a (z) 
is bounded inside the disk, since near a 



*(z) = 0((z - a) - *), F(z) = 0({z - a) - *). 



(4-55) 
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Thus the product x S(z)(F(z) © 1) may at most have square root singularities: however, comparing the 
RHPs that \& and F solve, we see that the product is a single- valued matrix, thus it must be analytic since 
at worst it may have singularities of type {z — a)~i . Being single-valued, however, it must have a Laurent 
series expansion that must , in fact, be a Taylor series by the previous a priori estimate of its growth. 
Therefore the product is actually analytic. Since F(z) is solely responsible for the unboundcdness of 1Z, 
this proves that V a {z) is bounded. 

4.6 Asymptotics of the biorthogonal polynomials 

Let us go back to modifications discussed in Section 4.1 of the RHP Problem 2.1 for Cauchy biorthogonal 
polynomials. It follows from (4-2), that within the half-lenses depicted in Fig. 4, 5 (on the ± sides of the 
cuts) the original matrix T{z) reads 



r±(aO 



1 

±e w i 1 



T ± {x) = C 1 T{x) ± G± 1 C- 1 



1 



1 



, x e A 



x€B. 



±e NV 2 1 

According to (2-9), the degree n monic biorthogonal polynomial is 

N 

r n ±e-T¥>i± ±r 12± eT^± 



(4-56) 
(4-57) 



Pn (x)=r n (x)=eH^+»^) 



-(v 1 -fl< 2 )+ 7+ ) 



{x e A). 



(4-58) 



Note that for x > 0, q( 2 \x) e M and has no jumps, and for x £ A we have ip x ± — — 2&#± {x) = 
±2in ( f* Q pi(s)ds - l) . From (4-2)-(4-9), 

T(z) = diag(l, -1, l)rp)diag(l, -1, 1) (zeC), T 12± {x) = ±T llT (x) = ±f 11± (x) (x e A), (4-59) 



and thus 



rue 



iNn /<Tn Pl( s ) ds 



p„( x)= e^( yi -B (2, +^)25R 
Vice versa, for z £ A (and outside of the right lenses) one has 

p„(z) = e A, s <1) r 11 ( 2 ) 



x £ A. 



(4-60) 



(4-61) 



Recall that T = (l + ©(iV^ 1 )) T , where T is defined in (4-11). Then one can obtain uniform asymptotic 
information on the behaviour of p n in any compact set of the complex plane. In particular, away from the 
endpoints, the asymptotics is expressible in terms of Theta functions (genus > 1) or algebraic expressions 
(genus 0). 

To obtain asymptotic information for the biorthogonal companions q n , one simply interchanges the 
roles of the measures. 
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5 Asymptotic spectral statistics and universality 



Once we have obtained a uniform asymptotic control of the biorthogonal polynomials, it is natural to 
investigate, in parallel to what has been done for the Hermitean matrix model (see, e.g. [10]), the large 
N behavior of the correlation functions for the Cauchy two-matrix model studied in [7] . We recall that 
the finite-size correlation functions for the spectra of M 1; M 2 distributed according to (1-1) correspond 
to a multi-level determinantal point process with the kernel 

3=0 3 

e -N(V 1 (x) + V 2 (-y)) 

Pj(x)Qi{y) 7 = hjSij , (5-2) 

2 x -f" y 

where {pj(x), Qj(y)}jeN are the monic biorthogonal polynomials. 2 
Define four auxiliary kernels 

H 00 {x,y) :=K N (x,y) , (5-3) 
r p -iWi(x)j T 

H 10 (y',y):= — - —H m {x,y) , (5-4) 

Jr+ y + x 

H m (x,x'):=l ( — H 00 (x,y) , (5-5) 

Jr+ y + x 

Remark 5.1 To obtain kernels for the correlation functions auxiliary kernels need to be multiplied by 
appropriate exponentials involving the potentials V\ , V% ■ 

Proposition 5.1 (Proposition 3.2 in [7]) 3 The auxiliary kernels are given in terms of the solution of 
the RHP in Proposition 2.1 as follows (for x,x' ,y,y' > 0) 

„ , , 1 [r-H-yMx)]^ x [v-\-y)T{-y%, 2 

H o(x,y) = - /n . s 2 , , H 10 {y ,y) = — , 

(z«7r) z x + y Zitt y y 

„ , n l [r-'VWh! „ , , p-^m-y)]^ 
H 01 (x 1 x) = — , Hii{y,x) 



(5-7) 



2i7r x — x' x + y 

One should not expect that the kernels Hqq(x, y), Hn(y, x), x, y > display any universal behavior, even 
in the scaling regime x = jf^ , y = The only "interaction" between the two matrices Mi and M2 
that can lead to a universality class is near the zero eigenvalue. 

On the other hand, if one considers separately statistics of the eigenvalues of Mi or M2, no new 
universality phenomena will appear, as we briefly explain below. It is sufficient to consider Hqi(x,x') 
since the computation for Hio(jj,y') is completely analogous. Again, we consider only the regular case. 

2 Normalizing factors hj appear in the definition of Kj\r(x, y) in contrast to the formula (3.31) in [7] since here we use 
monic rather than orthonormal biorthogonal polynomials. 
3 We adjusted the formulae from [7], see footnote 1. 
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Universality in the bulk. Let c belong to the interior of some cut. For simplicity we will consider 
only A, but the result extends to the other case with obvious modifications in the roles of the potentials. 
Due to regularity assumptions, at c we have pi(c) — C > 0. Consider 

X = C+ (?N> X ' = C+ (JN- M 
Then a straightforward computation using Proposition 5.1 yields 

pi{c)N 

( — { e %(<Pi+(*)-<Pi + (*')) _ e -f ! + {x) - V ! + {'■>:' ))\ (J + otN^ 1 )) . (5-9) 

2-7ri(£ — rj) V / 

Recall that 

/•CO 

(p 1 + (x) = i^sip 1 + (x) = — 2i7r / pi(s)ds (5-10) 

J x 



so that 



lim -JL- ^(x.z'Je-*^^^^ = e-&^ S ^ ~ f =: j^, ,) (5-11) 

AT->-oo p 1 [C)N 7r(£ — ?/) 

v v ' 

The prefactor above to the usual sine-kernel K S i n (^,r)) is not universal (it depends on the densities), 
however since K s [ n is conjugate to K sln , 

dct K sin (£i, = det /Wi',.;",.; . (5-12) 

Therefore all spectral statistics, gap probabilities etc. for Mi or M2 separately will follow the standard 
universality results for the Hermitean matrix model. 

Universality at the edge. Similarly, one finds near the edge x = a 

L ' n C:=lim^£, (5-13) 



CNi ' Cm ' ' «->a- v 7 ^ 

lim ^H 01 (x,x>)e-^W'» = ~ . (5 _ 14) 
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A An example: "double Laguerre" biorthogonal polynomials 



We consider the case where the two potentials are 
the same V\{x) = V2(— x) and are of the simplest 
possible form 



C -Vl(x) = x a e -bx 



(A.1) 



where both a, b > 0. 

We can rescale the axis and set b = 1 without 
loss of generality. 

The curve (3-4) appearing in Theorem 3.2 was 
computed in [2, 7] : 




y 



1 c 



2a 2 + 6a + 1 

3z 2 



27'=° 



According to section 4.3.1 we find a rational uni- 
formization of this curve as 

z=(l + a)t + 
2a + 1 



Figure 12: Depiction of the three sheets and how 
they are mapped onto three disjoint regions of the 
t-plane. The real cc-axis is mapped to the real t- 
axis and the two sides of each cut are mapped to 
(A. 2) the boundaries of the oval-shaped regions. The in- 
tersection of the ovals with the real i-axis are the 
four ramification points of the map z(t). 



V 



2a + 2 
1 



1 



1 



t - 1 



t + l 
1 



(A.3) 
(A.4) 



2a \(a + l)t-a (a + l)t + a / 

For a > there are four symmetric branch points on the real axis and the inner ones tend to zero as 
a 0, whereas all four tend to infinity as ±(a ± 2^fa) + 0(1) as a — > oo. Explicit formula? for the 



parametrix can be obtained by substituting u = 1 + a, U\ — u 2 

B Notation and main tools 



2a+l 
2a+2 



into (4-41). 



For a given smooth compact curve C of genus g with a fixed choice of symplectic homology basis of a 
and /3-cycles, we denote by uig the normalized basis of holomorphic differentials 



wit — 5jt , 

We will denote by G the Theta function 

9(z) := £ 

The Abel map u : C — >• C 9 with a base-point po is 

u(p) 



n-Tn—2t7rz-n 



\ r 




/ wi,.. 


> [ W 9 


-Jpo 





(B.l) 



(B.2) 



(B.3) 
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and is denned up to the period lattice Z + r • Z. For brevity we will omit any symbolic reference to the 
Abel map when it appears as argument of a Theta function: namely if p £ C is a point and it appears as 
an argument of a Theta-function, the Abel map will be implied, meaning that 

0(p — q) stands for 9(u(p) — u(q)) . 



We denote by K. the vector of Riemann constants 



9 V 



ue{p) 



' 33 



Pa 



(B.4) 



where the cycles oij are realized as loops with basepoint po and the inner integration is done along a 
path lying in the canonical dissection of the surface along the chosen representatives of the basis in the 
homology of the curve. 

The Riemann constants have the crucial property that for a nonspecial divisor V of degree g, T> = 
J2j=iPj> then tne function 

f(p) = Q(p-V-K) (B.5) 

has zeroes precisely and only at p = pj, j = 1 ... g. 

We also use Theta functions with (complex) characteristics: for any two complex vectors e, 5 the 
Theta function with half-characteristics e, S is defined via 



6 



, . e • t • e 1 _ 1 
(zj := cxp 2z7r ( — h -e • z 



2 4 

Here the half-characteristics of a point are defined by 

2z = 5 + re*. 

This modified Theta function has the following periodicity properties, for A,fieZ 9 



e 



(z + A + r/i) = cxp «7r(e*- A — S ■ p) — iirp ■ t ■ p — 2inz ■ p 



e 



(z) 



(B.6) 



(B.7) 



(B.8) 



Definition B.l The prime form E{p,q) is the (—1/2, —1/2) bi-differential on C x C 

®A(p-q) 



E(p,q) 



h.A(p)hA(q) 

9 

h A (p) 2 ln9 ^ 



fc=i 



uj k (p) =: uj A (p) , 



u=o 



where A = 



a choice of A. 



(B.9) 
(B.IO) 



is a half-integer odd characteristic (i.e. e ■ S is odd). The prime form does not depend on 
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The prime form E(p, q) is antisymmetric in the argument and it is a section of an appropriate line bundle, 
i.e. it is multiplicatively multivalued on C x C : 



E(p + aj,q) = E(p,q) , 



(B.11) 



E(p + Pj,q) = E(p,q)exp 




f 



) 



(B.12) 



In our notation for the half-characteristics, the vectors e, 5 appearing in the definition of the prime form 
are actually integer valued. We also note that the half order differential h A is in fact also multivalued 
according to 



Given a meromorphic function F (or a section of a line bundle) we will use the notation (F) for its divisor 
of zeroes/poles. For example, 



means that F has at most a pole of order k at p and a zero of multiplicity at least m at q. 
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